The role of quench rate in colloidal gels 
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Interactions between colloidal particles have hitherto usually been fixed by the suspension com- 
position. Recent experimental developments now enable the control of interactions in-situ. Here 
we use Brownian dynamics simulations to investigate the effect of controlling interactions on gela- 
tion, by "quenching" the system from an equilibrium fluid to a gel. We find that, contrary to the 
normal case of an instantaneous quench, where the local structure of the gel is highly disordered, 
controlled quenching results in a gel with a higher degree of local order. Under sufficiently slow 
quenching, local crystallisation is found, which is strongly enhanced when a monodisperse system is 
used. The higher the degree of local order, the smaller the mean squared displacement, indicating 
an enhancement of gel stability. 



I. INTRODUCTION 

Gelation is among the most striking features of soft 
matter [IH3]. Although many everyday materials are 
readily classified as gels, from toothpastes to yoghurts, a 
deep understanding of the gel state remains a challenge. 
In particular, three questions that might be tackled are, 
(i) what distinguishes a gel from a glass, (ii) which ma- 
terials can form gels, (iii) for a given material, what are 
the requirements for gelation? This work, which explores 
controlled quenching in colloidal gels with Brownian dy- 
namics simulations is relevant to (iii). 

For the first question, the identification of gelation with 
the crossing of a liquid-gas spinodal HJ [5] provides a 
working definition to distinguish gels from other dynam- 
ically arrested states, namely glasses. For the second, it 
seems that an effective attraction is required. Systems 
with short-ranged attractive interactions (up to around 
10% of the particle size) are well known to undergo gela- 
tion [1-3 . In the case of longer-ranged attractions in 
which the liquid state is stable, shallow quenches below 
the gas-liquid spinodal lead to complete phase separa- 
tion 0. Systems with short-ranged attractions and 
long-ranged repulsions can also exhibit similar behaviour 
I8UTU] . However the long-ranged repulsion significantly 
complicates the energy landscape, leading for example to 
degenerate mesophases [TTHT3] . Hereafter, we focus on 
systems without long-ranged repulsions, so the gels we 
consider are explicitly thermodynamically met ast able. 

Almost all gel literature concerns soft materials, that 
is to say, multicomponent systems of one or more meso- 
scopic component (colloids, nanoparticles, polymers) sus- 
pended in a fluid. However, gel-forming systems are often 
treated by integrating out the degrees of freedom of the 
smaller components (solvent molecules, small ions, poly- 
mers, etc) and considering an effective one-component 
system [14] . A reasonable question then is - can true one- 
component (molecular) systems form gels? Although to 
our knowledge no experiment has yet been performed, 
6*60 exhibits a property associated with a good gel- 



former, namely an attraction whose range is short com- 
pared to the molecular size. Cqq is predicted to form a 
gel [15] , and even the longer-ranged Lennard- Jones model 
(relative to the molecular diameter) will form a gel under 
a sufficiently deep quench |B]. 

This brings us to the third question - the requirements 
for gelation in a given material which forms the subject 
of this article. The gels we consider are intrinsically non- 
equilibrium, and therefore, the means of preparation is 
important. Here we shall consider a system with a rel- 
atively short-ranged attraction, in which the equlibrium 
state is gas-crystal phase coexistence. Forming a gel, 
therefore, means avoiding phase separation to a gas and 
a crystal. Moreover, the resulting state must be dynami- 
cally arrested, which also requires that it percolates [16J. 
Possible routes to gelation are illustrated in Fig. [T] The 
most obvious route is quenching, reducing the (effective) 
temperature. Quenching must be carried out with suf- 
ficient speed that phase separation does not occur. For 
soft matter systems, this is often straightforward. 

In colloids for example, where the absolute tempera- 
ture is typically held at 298 K, effective temperature is 
varied by changing the interaction strength. The effective 
temperature is then related to the depth of the attractive 
well of the interaction potential, and is often fixed for a 
given sample. Scanning a phase diagram then requires 
preparation of a considerable number of different sam- 
ples, each with its own effective temperature. Following 
preparation of a sample, for example by shear, the sys- 
tem may be said to have undergone an (uncontrolled) 
instantaneous quench. A typical example is a colloid- 
polymer mixture, where the effective temperature is set 
by the polymer concentration (through the depletion in- 
teraction) [17]. In this case, the interactions are mediated 
by the polymer on much shorter timescales that the larger 
(and slower) colloids, and it is reasonable to consider an 
instantaneous quench as a route to gelation. 

At smaller lengthscales, in nanoparticle and molecular 
systems, two main effects come into play. Firstly interac- 
tions are often reasonably constant over the temperature 



range of interest, but temperature is used as a control pa- 
rameter. A consequence then is that quench rates cannot 
be instantaneous is often the case for colloids. This is sig- 
nificant, as in the case of small nanoparticles and molec- 
ular systems, the dynamics are often fast enough that 
quench rates which avoid phase separation can be tech- 
nically challenging [15]. In this case, crunching, where a 
gas of clusters is compressed, remains a route by which 
gels can be realised. 

However, since gels are usually formed of soft matter 
with mescoscopic components, it is natural that most 
work on gelation assumes that the interactions between 
the particles were fixed. In other words, that quenches 
were instantaneous. This is typically the case in both 
experiments p] [3] and computer simulation [TJ [TDJ H8l - 
[20] . Recent developments in controlling colloid-colloid 
interactions allow the interactions (and thus the effective 
temperature) to be changed at will, even on timescales 
much faster than the colloid dynamics, such that control- 
lable quenches may be carried out [21]. It has also be- 
come possible to control attractive interactions between 
colloids such as temperature-dependent depletion attrac- 
tions [22j[23], multiaxial electric fields [24] and the critical 
Casimir effect [25H27] . This opens the way to consider 
the role of controlling effective temperature (by chang- 
ing interactions) in an experiment and thus enables us 
to consider the role of quench rate in colloidal gelation. 
Varying the quench rate has been carried out in gela- 
tion in molecular systems [15 . There, using molecular 
dynamics simulations, a high quench rate was found to 
be necessary to prevent phase separation to equilibrium 
gas-crystal coexistence. Here we shall use Brownian dy- 
namics simulations to model colloidal gelation in which 
we shall consider the effect of quench rate upon the gels 
formed. 

This paper is organised as follows. First we intro- 
duce the model system, which is a representation of a 
colloid-polymer mixture using Brownian dynamics simu- 
lations [28] . We then consider the response of the system 
to a conventional treatment of an instantaneous quench 
where the effective temperature is fixed at the outset. 
The system is characterised through a novel analysis of 
local structure we have developed, the topological clus- 
ter classification (TCC) [23 ESI. The effect of quench 
rate is then presented in both polydisperse and monodis- 
perse systems. Finally we discuss the implications of our 
findings. 



II. METHODS 

A. Model 

We model a colloid-polymer mixture with Brownian 
dynamics simulations. In the experimental system upon 
which we base our simulations, residual electrostatic in- 
teractions are screened [28, 30 and are neglected here. 
We found good agreement with simulations using the 
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FIG. 1: (colour online) Routes to gelation in short- 
ranged attractive systems. Gelation results upon quench- 
ing, provided the quench rate is sufficient to avoid phase sep- 
aration (shaded region). Crunching of isolated clusters pro- 
vides an addition route to form a gel where accessible quench 
rates are too slow to prevent demixing. Here "C" denotes 
isolated clusters, "G" gas, "F" (supercritical) fluid. 



Morse potential [28 . The Morse potential reads 



f3u(r) = Pee- p0 ^- r) ( e -P^~ T ) - 2) 



(1) 



where j3 = l//c#T the thermal energy. The potential is 
truncated and shifted at r = 1.4a. The effective tem- 
perature is controlled by varying the well depth of the 
potential e. We set the range parameter po = 25.0. Col- 
loidal particles are typically polydisperse, and here we 
treat polydispersity by scaling r in Eq. [I] by a Gaussian 
distribution in a with 4% standard deviation (the same 
value as the size polydispersity in the experimental sys- 
tem). In addition we consider the case of a monodisperse 
system. We express time in units of the Brownian time 
i~b = ( cr /2) 2 /6D which is the time taken for a particle 
to diffuse its own radius. Here D is the diffusion con- 
stant. In the simulations, tb ~ 894 time units. The 
time step is 0.03 simulation time units. Where the ef- 
fective t emper ature is fixed throughout the simulation 
(Section |III A), the runs are equilibrated for 5.0 x 10 6 



steps and run for a further 5.0 x 10 b steps. These sim- 
ulation runs therefore correspond to approximately 168 
Brownian times or around 10 min, a timescale certainly 
comparable to experimental work. Except in the case of 
long quench times (1.01 x 10 4 tb) where 2048 particles 
were used, the system size was 10000 particles. We fix 
the packing fraction <j) = 7rpcr 3 /6 = 0.1 throughout. At 
this packing fraction, quenching below the critical point, 
which is approximately located at (f3e c « 2.69) [31 , leads 
to a percolating network - i.e. a gel. The critical point 
gives a reasonable estimate for gelation, since the spin- 
odal line is found to be very flat for short-ranged systems 
[32] . In the case of simulations where a quench is carried 
out, j3e is increased linearly from unity to a chosen value, 
at which point the system is run for 5 x 10 6 time steps 
and analysed. 
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FIG. 2: Clusters detected by the topological cluster 
classification. These structures are minimum energy clus- 
ters of the Morse potential with p = 25.0 [33] . 



B. Structural Analysis - Topological Cluster 
Classification 

To analyse the structure, we identify the bond network 
using the Voronoi construction with a maximum bond 
length of 1.4 a. Having identified the bond network, we 
use the Topological Cluster Classification (TCC) to de- 
termine the nature of the cluster [29, 34 i . This analysis 
identifies all the shortest path three, four and five mem- 
bered rings in the bond network. We use the TCC to find 
clusters which are global energy minima of the Morse po- 
tential for the range we consider (po = 25.0), as listed in 
[33] and illustrated in Fig. [2j In addition we identify the 
thirteen particle structures which correspond to FCC and 
HCP in terms of a central particle and its twelve nearest 
neighbours. For more details see [29]l34]. 



III. RESULTS 

The results section is organised as follows. We begin 
by presenting a brief overview of the effect of quench- 
ing. We consider the conventional case of instantaneous 
quenching, and provide an analysis from the perspective 
of local structure. We next investigate a finite quench 
rate, and the effect of varying quench rate on the result- 
ing structure and dynamics of the gel. Having established 
the general nature of quenching on the local structure, 
we examine the special case of a monodisperse system. 
Dynamic data is then employed to investigate relative 
stabilities of the two routes to gelation - instantaneous 
quenching and quenching at finite rates. 

An overview of the system is shown in Fig. [3] These 
snapshots show the network structure that is developed 
by colloidal gels at a packing fraction <\> = 0.1 and at- 
traction j3e = 5.0. While both are networks, the detail 
of the local structure is markedly different, as revealed 



by the topological cluster classification. The instanta- 
neously quenched (polydisperse) system [Fig. [3^a)] is 
dominated by a 10 membered fivefold symmetric local 
structure, while the monodisperse system with a quench 
time 671 tb shows a significant degree of local crystallisa- 
tion but still retains a similar overall network structure. 
A slight degree of coarsening is apparent in the static 
structure factor S(k) in Fig. [3Fc) relative to the instan- 
taneous quench. We also note that, for the same quench 
time, a monodisperse system appears to coarsen some- 
what more than a polydisperse system. However, the 
change in S(k) appears much less dramatic than that in 
local structure illustrated in Fig. |3ja) and (b). 



A. Instantaneous quenches in polydisperse systems 
- a local structural analysis 

We now examine instantaneous quenches in more de- 
tail. A conventional structural analysis is presented in 
Fig. pFa) in the form of the static structure factor S(k). 
This illustrates the typical behaviour of systems undergo- 
ing dynamical arrest: the structure in the fluid ((3e < 3) 
responds to changes in effective temperature, in that the 
low k region of S(k) increases strongly as phase separa- 
tion is approached. Conversely, upon dynamical arrest 
for j3e > 5, no further change in the static structure is 
seen. Thus one might conclude that under instantaneous 
quenching, all gels with <\> = 0.1 have the same local 
structure. 

In fact this is not so, as shown in Fig. H[b). Here the 
topological cluster classification is shown for the same 
data. This analysis reveals the generic nature of gela- 
tion at the single-particle level: at high effective tem- 
peratures, in the ergodic fluid phase, few clusters are 
found. Those few that are present are predominantly 
five-membered triangular bipyramids. Upon gelation, 
the cluster population rises very sharply. This indicates 
that the particles condense into clusters which comprise 
the gel network [Fig. J3^a) and (b)]. However it is the 
type of the clusters which reveal the local structure of 
the gel. Although the equilibrium state is gas-crystal co- 
existence, even at the level of the few particles associated 
with a cluster, the kinetic pathway to crystallisation is 
arrested and very few particles ae found in a local HCP or 
FCC environment. Upon deeper quenching, the system 
moves further from equilbrium, and even fewer particles 
in local crystalline environments are found. Rather than 
crystalline environments, particles are found in clusters of 
between 8 and 10 particles. As illustrated in Fig. [2j these 
are built around pentagonal bipyramids and are thus five- 
fold symmetric. Such five-fold symmetry has long been 
believed to suppress crystallisation [35]. We note that 
these same clusters are prevalent in dense fluids of both 
Morse particles and hard spheres [55] , This feature of 
crystallisation suppressed by five-fold symmetric clusters 
is found also in experiments on colloidal gels [28], simu- 
lations of molecular gels [15], and isolated clusters [37] . 
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FIG. 3: Simulation snapshots labelled by the topological cluster classification, (a) Polydisperse system instanta- 
neously quenched, (b) Monodisperse system with a quench time of 671 tb- Both (a) and (b) correspond to an identical state 
point with /3s = 5.0. Color of particles on right of (b) denotes cluster in which the particle is identified by the TCC as shown 
in Fig. [2] Note the vastly increased population of particles in crystalline environments in the case of the monodisperse system 
with a quench time of 671 tb- (c) Static structure factor S(k) for state points shown in (a) (purple line) and (b) (pink line) 
along with a polydisperse system with a quench time of 671 tb (navy line). 
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FIG. 4: The conventional situation: instantaneous quench, (a) The static structure factor (S(k)) for varying strengths 
of attraction. Where gelation occurs on the simulation timescale (f3s > 3), the S(k) are indistinguishable for practical purposes. 
(b) The topological cluster classification reveals changes in local structure upon quenching. Dashed line denotes the strength of 
attraction at criticality e c which approximately locates gelation. N c is the number of particles in a given cluster, N is the total 
number of particles, (c) mean squared displacement for the well depth f3e (inverse effective temperature) indicated. Dashed 
grey line has a slope of unity (diffusive motion). 



where in the case of the latter, the ground state is one 
of the clusters considered by the TCC. However, sup- 
pression of crystallisation is not universal, experiments 
in which locally crystalline clusters are found have also 
been reported [38] . 

The change in dynamics upon gelation for instanta- 
neous quenches is shown in the mean squared displace- 
ment (MSD) data in Fig. pFc). Two observations are 
important concerning this plot. The first is that all state 
points are diffusive. There is little indication of a plateau 
that one might associate with slow dynamics. The second 
is that although in the fluid state, the system is rather 
insensitive to the effective temperature with only a slight 
decrease in MSD upon increasing the strength of the at- 
traction. Upon gelation there is a sudden drop in mobil- 
ity of more than an order of magnitude. Although MSD 
indicates diffusive motion, the magnitude of the motion is 
small, only reaching the particle lengthscale after around 
IOOtb- As significant as this drop in mobility is, some 
restructuring cannot be ruled out, indicating that these 
gels age. 



B. Finite quenches in polydisperse systems - 
structure 



Let us consider finite quench rates. We shall use our 
local structural analysis, the TCC. As might be sur- 
mised from Fig. [4^a), the static structure factor is not 
very sensitive to quench rate. TCC analyses of data for 
quench times of 671 tb and 1.01 x 10 4 t# are shown in Fig. 
[5Va) and (b) respectively. In the corresponding experi- 
mental system, the quench times are around 30 minutes 
and 8 hours respectively. For comparison, instantaneous 
quench data is plotted in Fig. [5] (c). Here we quench 
from a fluid at f3e = 1 to the final value of j3e as shown 
in Fig. [5|a) and (b). Each effective temperature there- 
fore corresponds to a different simulation run. We begin 
by considering the moderate quench time of 671t#. The 
generic behaviour is similar to that of an instantaneous 
quench [Fig. (5^c)], with a sharp rise in cluster population 
upon gelation (the very small cluster population in the 
equilibrium fluid (3e < j3e c is of course unaffected by the 
quench protocol). For shallow quenches (/3e ~ 5), there is 
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FIG. 5: The effect of a quench rates on the local structure in polydisperse systems (a) TCC analysis for a quench 
time of 671Tb- (b) TCC analysis for a quench time of 1.01 x 10 4 tb. (c) Data from an instantaneous quench. Dashed lines 
denote e c 



a qualitative similarity: the structure is dominated by the 
10-membered cluster. However, upon quenching deeper a 
significant difference emerges. The effect is considerable 
and is summarised as follows: In the case of a moderate 
quench time, the fivefold symmetric 10-membered cluster 
population shows a rapid rise upon gelation up to a pop- 
ulation of around N c /N « 0.42 and continues to increase 
slowly upon deeper quenching, finally reaching a value of 
N c /N « 0.5 for f3e = 20. Very different behaviour is seen 
in the case of the instantaneous quench. Although the 
10-membered cluster again exhibits strong rise upon gela- 
tion, its population falls upon deeper quenching, finally 
reaching a value around N c /N « 0.3 for (3e = 20. Other 
smaller clusters, notably eight- and nine-membered clus- 
ters form a very significant fraction of the population in 
the case of deep quenches. 

For a long quench time (1.01 x 10 4 t#), Fig. [H|b) 
shows a further change in local structure. Higher-order 
clusters again are favoured. In particular, for 4 < j3e < 5, 
the most popular cluster is the HCP crystalline environ- 
ment, with a value of N c /N ^0.3. In other words, 
around 30% of the system has succeeded in reaching, 
or at least getting close to, local equilibrium. Even 
at deeper quenches, although the crystal population is 
smaller (N c /N « 0.15), much more of the system has 
crystallised compared to shorter quench times. The rise 
in HCP environments is accompanied by a rise in the 
13-membered five- fold symmetric cluster. This pentag- 
onal prism which is the minimum energy structure for 
13 Morse particles is similar to - but distinct from - the 
icosahedron. 



Finite quenches in monodisperse systems - 
structure 



Polydispersity is known to suppress crystallisation. We 
therefore consider the effect of setting the polydispersity 
to zero (a monodisperse system), keeping all other condi- 
tions the same. The effect on the local structure is shown 
in Fig. |6j Except that we now consider a monodisperse 
system, the situation is identical to Fig. [5] Perhaps the 



most striking result is the instantaneous quench, shown 
in Fig. |6Fc). The similarity to the polydisperse case 
is remarkable, and crystallisation is strongly suppressed. 
This is consistent with the idea that under instantaneous 
quenching, very little re- arrangement is possible, parti- 
cles remain kinetically trapped in the clusters they form 
upon condensation, and lack the thermal energy to rear- 
range. 

The case of the finite quenches [671 tb and 1.01 x 10 4 
tb in Fig. [6Fa) and (b) respectively] is very different in- 
deed. Unlike the dominance of the five-fold symmetric 
10-membered cluster in the polydisperse case, here the 
crystalline HCP and FCC clusters are the most abun- 
dant. Like the polydisperse system, the degree of crys- 
tallisation appears to be increased in the case of slow 
quenching (1.01 x 10 4 r#). As shown in Fig. |3[b), this 
monodisperse system can largely reach local equilibrium, 
however it retains a network structure: a crystalline gel. 



D. Finite quenches - dynamics 

The effect of quench times upon the dynamcs is shown 
in Fig. [7| As in the instantaneous case, the mean squared 
displacement appears diffusive, with no indication of a 
plateau. However, as noted above, the magnitude of the 
MSD is small, with the exception of j3e = 3. For the 
finite quench times, the MSD drops significantly com- 
pared to the instantaneous quench time. For the long 
quench time, (1.01 x 10 4 r#), this drop is very consider- 
able indeed, with a reduction of more than an order of 
magnitude. At the very least, this suggests that the finite 
quench time leads to a gel with a much higher stability 
Similar effects are found in the case of monodispersity, 
which seems consistent with the increase in crystalline 
order observed above. 



IV. DISCUSSION 

We have shown that colloidal gels quenched slowly ex- 
hibit local structural differences compared to the same 




FIG. 6: The effect of a quench rates on the local structure in monodisperse systems (a) TCC analysis for a quench 
time of 671Tb- (b) TCC analysis for a quench time of 1.01 x 10 4 tb. (c) Data from an instantaneous quench. Dashed lines 
denote e c 




FIG. 7: The effect of a quench times on the dynam- 
ics, (a) Polydisperse system (b) monodisperse system. Mean 
squared displacement for quench depth and quench times. 
Solid lines - instantaneous quench, dotted lines, moderate 
(671tb) quench time, dashed lines, long (1.01 x 10 4 tb) quench 
time. 



system quenched instantaneously (in other words, where 
the particle interactions are fixed). The local structural 
motifs of instantaneously quenched gels are are five-fold 
symmetric, and are also found in dense equilibrium fluids 
for this system [36]. It seems reasonable to suppose that, 
upon quenching, the system begins to phase separate, 
and the dense fluid phase - whose arrest drives gelation 
- is structurally similar to high-density fluids. In partic- 
ular, the local structure is dominated by fivefold sym- 
metric clusters whose principle motif is a five-membered 
ring (7 < m < 10 in Fig. |2l). Perhaps surprisingly, the 
local structure upon instantaneous quenching is rather 
insensitive to polydispersity 



A. Structure and finite quench rates 

Quenching at a finite rate on the other hand gives the 
system some chance to reach thermodynamic equilibrium 
as it is "cooled" . Specifically, in equilibrium, phase sepa- 
ration is expected for e > e c . As mentioned in the intro- 
duction, here we are concerned with a system with short 
range attractive interactions in which the gel is explicitly 



metastable to crystal-gas coexistence. However, even for 
a quench rate of 4.0 x 10 -4 Per^ 1 , no such phase separa- 
tion is observed, although some coarsening is observed for 
a moderate quench time of 671 tb- In the case of exper- 
iments on colloid-polymer mixtures, due to the change 
in polymer size with temperature 1°C corresponds to a 
change in attraction of order /c#T. The resulting change 
in effective temperature (due to the polymer-induced de- 
pletion attraction) corresponds to a quench rate of just 
~ 10 -3 °Cs _1 . Thus the overall feature of colloidal gela- 
tion is remarkably robust to temperature quenching. Or, 
equivalently, colloidal model systems are very slow to 
equilibrate. 

Although the system continues to gel, even with mod- 
erate quench rates (quench time 671 7b), the local struc- 
ture is a strong function of cooling rate. In particular, 
the structure of the gel is rather constant with quench 
depth [Fig. |5ja)]. In other words, in the early stages of 
gelation, soon after the system has become metastable 
and begin to phase separate, the particles assemble into 
10-membered clusters. Upon deeper quenching, this lo- 
cal structure is preserved, as the particles gradually lose 
their mobility. By contrast, systems subjected to in- 
stantaneous quenches are less able to self-organise into 
such large clusters and, particularly at deeper quenches, 
a substantial number of smaller clusters are present. Of 
particular note is the 5-membered triangular bipyramid. 
This is formed from two tetrahedra, and the tetrahedron 
is the simplest rigid assembly of spheres in 3D. In other 
words, instantaneous quenches lead to very simple struc- 
tures, which suggest rapid formation with no further re- 
arrangement, such as is the case for diffusion-limited clus- 
ter aggregation. 

The opposite case then, is a slow quench (1.01 x 10 4 
tb)- Here, although full phase separation does not oc- 
cur, a very considerable number of particles are able to 
reach local equilibrium and crystallise. We have based 
our discussion on a slightly polydisperse system as found 
in experiments [28, 30 . This has two main consequences. 
Firstly, crystallisation should be suppressed. Although 
hard spheres with a polydispersity of 4% as we have em- 
ployed here do crystallise, short-ranged attractive sys- 



terns are known to be more sensitive to polydispersity 
[2]. Secondly, the relationship between polydispersity 
and the eventual approach to equilibrium is complex. 
What is known is that polydisperse hard spheres phase 
separate into less polydisperse daughter populations [39] , 
and it seems reasonable to suppose that similar behaviour 
might occur here. At the very least, this suggests that full 
equilibration would take a long time. Conversely, in the 
monodisperse case, a far higher degree of crystallisation is 
found, even though for instantaneous quenches the local 
structure seems insensitive to polydispersity. However we 
note that previous work did indeed find some crystallisa- 
tion in Brownian dynamics simulations of monodisperse 
attractive particles where quenching was instantaneous, 
but this tended to be supressed by gelation - crystallisa- 
tion occured at weaker attractions than e c [20 . Such 
observations are consistent with measures of the self- 
assembly from a perspective of local reversibility as de- 
termined by deviations from the fluctuation-dissipation 
theorem [40J, where optimal assembly (crystallisation) 
occurs at e < e c . 



B. Dynamics and finite quench rates 



C. The role of dynamics 



Before closing, a few words on the role of the dynam- 
ics employed are in order. In short, this is the observa- 
tion that, while here we have used Brownian dynamics 
and seen no phase separation, in a comparable system 
with molecular dynamics (MD) [15], it was hard to pre- 
vent complete phase separation. While the systems were 
not accurately mapped to one another, (the MD system 
has a small region where the liquid is thermodynami- 
cally stable), the extent of the difference in behaviour 
seems to warrant further investigation, although some 
previous work did not reveal a very large difference be- 
tween the two forms of dynamics [41]. Can a system with 
Newtonian dynamics "know" its way through the energy 
landscape so much better than a system with Langevin 
dynamics that it is hard to kinetically trap the former 
and apparently impossible to avoid kinetic trapping in 
the latter? We hope to answer this question in the near 
future, and also to develop the results we have presented 
here with larger scale simulations that can explore the 
role of finite size effects. Finally, we note that solvent- 
mediated hydrodynamic interactions have been shown to 
have a profound effect upon gelation and may also influ- 
ence local structure [42] . 



The drop in MSD induced in the same system by vary- 
ing the quench rate is remarkable, and surprising in its 
extent. At one level, it underlines the need to use finer 
probes of structure than S(k) in order to unravel the 
origins of the dynamical behaviour. We can rationalise 
the reduction in MSD as consistent with the system being 
better equilibrated, i.e. lower in its energy landscape. As 
intuitive as this observation is, backed up by structural 
observations consistent with improved equilibration, the 
strong dependence of the dynamics is nonetheless worthy 
of further investigation. Moreover, the MSD data raises 
additional questions. 

For example, why are the MSDs diffusive? Although 
the sampling times are quite short, the absence of any 
plateau can be interpreted as indicative of further ir- 
reversible structural rearrangement. However, that all 
the MSDs have such similar slopes almost regardless of 
quench depth and time remains a curiousity. 

A further comment on the approach to equilibrium and 
stability concerns aging. While this has already recieved 
attention for instantaneuous quenching [19] [41], an im- 
portant question concerns the link between aging and 
quench rate. Certainly, instantaneously quenched sys- 
tems are higher in the free energy landscape and exhibit 
more mobility. Longer quench times produce systems 
lower in the energy landscape with less mobility, i.e. pre- 
cisely the consequences of ageing |3]. Can quenching be 
related to aging? Is it direct lty equivalent? Connecting 
our measures of local structure with dynamical quantities 
such as MSD provides a means to tackle these questions. 



V. CONCLUSIONS 



We have carried out Brownian dynamics computer sim- 
ulations of a colloidal gel. By "cooling" the system from 
the stable fluid to a metastable gel, we have explored the 
role of quench rate in colloidal gels. Our results indi- 
cate that, despite its metastable nature, gelation is sur- 
prisingly robust. All quench rates we have been able to 
perform resulted in a gel, none exhibited phase separa- 
tion. Conversely, the local structure accessed through 
the topological cluster classification is a strong function 
of quench rate. Decreasing the quench rate enables the 
system to access larger structures associated with lower 
basins in the energy landscape, and, for the slowest rates 
we access, a considerable part of the system, can reach 
local equilibrium and crystallise. Monodisperse systems 
exhibit a much higher degree of crystallisation than only 
slightly polydisperse systems. This is reflected in the dy- 
namics: the mean squared displacement drops markedly 
for the same state point when a slow quench is employed. 
We therefore conclude that local structure, as measured 
by the TCC, is strongly coupled to the dynamics. Our 
results suggest that controlled quenching, as is beginning 
to become possible in experiments on colloidal model sys- 
tems, may be a means by which colloidal gels can be stab- 
lised, and products based on gels may enjoy an extended 
shelf life. 
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